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Abstract 

We present a method for proving the existence of symmetric periodic, 
heteroclinic or homoclinic orbits in dynamical systems with the reversing 
symmetry. As an application we show that the Planar Restricted Cir- 
cular Three Body Problem (PCR3BP) corresponding to the Sun-Jupiter- 
Oterma system possesses an infinite number of symmetric periodic orbits 
and homoclinic orbits to the Lyapunov orbits. Moreover, we show the 
existence of symbolic dynamics on six symbols for PCR3BP and the pos- 
sibility of resonance transitions of the comet. This extends earlier results 
by Wilczak and Zgliczynski 12 . 

1 Introduction. 

The Planar Restricted Circular Three Body Problem (PCR3BP) has attracted 
much attention of scientists. In particular, some transport properties of this 
system may be applied in space mission design (see 0] and references given 
there). The problem has been studied by Koon, Lo, Marsden and Roos in 
4 , where the numerical evidence of the resonance transitions for the PCR3BP 
for the parameter values corresponding to the Sun-Jupiter-Oterma system is 
presented. The rigorous proof of some facts discovered in ^ was given by 
Wilczak and Zgliczynski in |^] (see also Stoffer and Kirchgraber paper 
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In the present paper, as in J2]> we restrict our attention to the following 
parameter values for PCR3BP C — 3.03, /i = 0.0009537 - the parameter values 
for Oterma comet in the Sun- Jupiter system (see 0]). For this parameter values 
there are two hyperbolic periodic orbits L\ and L\, the Liapunov orbits, around 
the libration points L\ and L2, respectively. In |12| and [8] it was proven that 
there exist homoclinic solutions to both LI and L* 2 periodic orbits and a pair of 
heteroclinic connections between them in both directions. 

In this paper we present the proof of the following facts: 

• The PCR3BP possesses two pairs of homoclinic orbits both to L\ and L\. 
These homoclinic orbits are geometrically different. Informally speaking 
they are close to different resonances, namely 3:2, 5:3 for the orbits ho- 
moclinic to L\ and 1:2, 2:3 for the orbits homoclinic to L^. Moreover, it 
is possible for a comet to move between these four resonances in arbitrary 
order. 

• The PCR3BP possesses an infinite number of geometrically different sym- 
metric periodic and homoclinic orbits. 

Let describe now what constitute a new element in the present paper. While 
the numerical evidence of three of the above mentioned homoclinic orbits is given 
in 0], the 2:3 homoclinic orbit appear to be a new one. The technique of the 
proof of the existence of these homoclinic orbits and the symbolic dynamics is 
the same as in i- e - combines topological tools (covering relations) with 
rigorous numerics. 

The main novelty of the present paper, when compared to |^1 112j . is the 
existence of an infinite number of geometrically different symmetric periodic 
and homoclinic orbits. While the numerical evidence of the simplest symmetric 
homoclinic orbits is given in 0] the numerical method used there cannot yield 
the existence of an infinite number of them even with the help of validated 
numerics. In this paper we use the topological method introduced recently by 
Wilczak in and developed later by both authors in [E]. For the purpose of 
this introduction we briefly describe the main points of method for symmetric 
periodic points. The method is based on two observations: 

• to detect symmetric periodic orbits for map P (a Poincare map) with an 
reversing symmetry R (in the PCR3BP case composition of a suitable 
reflection and the time inversion) it is enough to look for intersections of 
Fix(i?) = {x I R(x) = x} with P k (Fix(R)). This is the Fixed Set Iteration 
method [HE] (also known as DeVogelaere method 3 ). Any point from 
such intersection give rise to 2&-periodic point. 

• covering relations give some control of pieces of P fe (Fix(i?)), which make 
it possible to prove that P fe (Fix(S'))nFix(i?) is nonempty for k sufficiently 
large and that the period of the periodic point is indeed equal to 2k. 

The paper is organized as follows. In Section[2]we recall the PCR3BP and its 
properties. In Section the proof of the existence of a new pair of homoclinic 
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orbits both in exterior and interior regions to the Lyapunov orbits L*, L\ is 
presented. In Section 0] the symbolic dynamics on six symbols is established. 
We also discuss the resonance transitions there. In Section [5] the existence of 
symmetric periodic and homoclinic orbits is proven. 

Throughout the paper we will use the definitions and notations from |12|. 



2 Short description of the system. 

We follow papers an( A use the notation introduced there. 

Let S and J be two bodies called Sun and Jupiter, of masses m s = 1 — /i 
and rrij — fj,, fj, € (0, 1), respectively. They rotate in the plane in circles counter 
clockwise about their common center and with angular velocity normalized as 
one. Choose a rotating coordinate system, so that origin is at the center of 
mass and the Sun and Jupiter are fixed on the x-axis at (— /x, 0) and (1 — fi, 0) 
respectively. In this coordinate frame the equations of motion of a massless 
particle called the comet or the spacecraft under the gravitational action of Sun 
and Jupiter are (see j3] and references given there) 

x - 2y = Cl x (x,y), y + 2x = Q y (x, y), (1) 

where 

x 2 + y 2 1-n fi fi(l - n) 



tt(x,y) 



n r 2 



r\ = \J (x + /i) 2 + y 2 , r 2 = y/ (x - 1 + fi) 2 + y 2 

Equations are called the equations of the planar circular restricted three- 
body problem (PCR3BP). They have a first integral called the Jacobi integral, 
which is given by 

C{x 1 y,x 1 y) = ~{x 2 + y 2 ) + 2n{x 1 y). (2) 

We consider PCR3BP on the hypersurface 

M(fi,C) = {(x,y,x,y) \ C(x,y,x,y) = C}, 

and we restrict our attention to the following parameter values C = 3.03, fi — 
0.0009537 - the parameter values for Oterma comet in the Sun- Jupiter system 
(see H). 

The projection of M([i,C) onto position space is called a Hill's region and 
gives the region in the (a;, y)-plane, where the comet is free to move. The Hill's 
region for the parameter considered in this paper is shown on Figure^in white, 
the forbidden region is dark. The Hill's region consists of three regions: an 
interior (Sun) region, an exterior region and Jupiter region. 

As was mentioned in the Introduction we restrict our attention to the fol- 
lowing parameter values C = 3.03, /i = 0.0009537 - the parameter values for 
Oterma comet in the Sun- Jupiter system (see 0]). Since we work with fixed 
parameter values we usually drop the dependence of various objects defined 
throughout the paper on n and C, so for example M. = M(fj,,C). 
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2.1 Poincare maps. 



We consider Poincare sections: 9 = {(x, y,x,y) € M \ y = 0}, 6 + = 6fl{j/> 

o}, e_ = en{y<o}. 

On 6± we can express y in terms of x and x as follows 
y = ±y/2fl(x,0) -x 2 - C 

Hence the sections 0± can be parameterized by two coordinates (a;, x) and we 
will use this identification throughout the paper. More formally, we have the 
transformation T± : M 2 — > 9± given by the following formula 

T±(x,x) = {x,0,x,±y/2n{x,0) - x 2 - C ) 

The domain of T± is given by an inequality 2fl(x, 0) — x 2 — C > 0. 

Let iT± '■ B± — > M. and n x '■ 0± — > K denote the projection onto x and a; 
coordinate, respectively. We have ir±(x ,x ) = x and ir x (x ,x ) = x - 

We will say that (x, x) € &± meaning that (a;, x) represents two-dimensional 
coordinates of a point on 6±. Analogously we give a meaning to the statement 
M C 9± for a set M C K 2 . 

We define the following Poincare maps between sections 



p+ 


:6+- 






: e_ - 


-» e_ 


2 


:6+- 


-» e_ 


Pi : 

2 ' 


e_ - 


>e+. 



As a rule the sign + or — tells that the domain of the maps P± or P* t ± is 
contained in 6± (the same sign). Observe that 

P+(x) = P h ._ o Pi i+ (aO, P-{x) = P h+ o P h _(x) 

whenever P+(x) and P-{x) are defined. These identities express the following 
simple fact: to return to 0+ we need to cross O with negative y (this is Pi + 
first and then we return to O with y > (this is Pi _). 

Sometimes we will drop signs in P± and Pi_±, hence P(z) = P+(z) if z e 0+ 
and P(z) = P-(z) if z e 6_, a similar convention will be applied to Pi. 

2.2 Symmetry properties of PCR3BP 

Notice that PCR3BP has the following symmetry 

R{x, y, x, y, t) = (x, -y, -x, y, -t), 

which expresses the following fact, if (x(t),y(t)) is a trajectory for PCR3BP, 
then (x(— t), — y(— t)) is also a trajectory for PCR3BP. From this it follows 
immediately that 

if P ± (.t ,.to) = (xi,±i) then P±(x 1 ,-x 1 ) = (x ,-x ), 
if Pi i± (x ,xo) = (xi,x\) then Pi- f (xi,-xi) = (x ,-x ). 
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We will denote also by R the map R : 9± — > 0± x) = (x, —x) for 
(x, x) € 8±. Now Eq. 10 can be written as 



3 The existence of new homoclinic orbits. 

The goal of this section is to present the proof of the existence of new homoclinic 
orbits with different resonances. 

The notion of the resonance. We rewrite here an informal definition of 
the resonance from Sec. 5.1]. Recall that the PCR3BP is a perturbation 
of the two-body problem. Hence, outside a small neighborhood of Jupiter, the 
trajectory of a comet follows essentially a two-body orbit around the Sun. In 
the heliocentric inertial frame, the orbit is nearly elliptical. The mean motion 
resonance of the comet with respect to Jupiter is equal to a -3 / 2 where a is the 
semi-major axis of this elliptical orbit. Recall that the Sun-Jupiter distance is 
normalized to be 1 in the PCR3BP. The comet is said to be in p : q resonance 
with Jupiter if a~ 3 / 2 w p/q, where p and q are small integers. In heliocentric 
inertial frame, the comet makes roughly p revolutions around the Sun in q 
Jupiter periods. Observe that this definition of the resonance make also sense 
for the orbits, which are non-periodic (for example orbits homoclinic to L\ or 
^2) j we J us t have to compute the semi- major axis for the piece of orbit away 
from Jupiter. A heuristic approach, which allows to read the resonance of an 
orbit from the trajectory in the rotating frame is described in Appendix. 

In |12| the following theorem was proved. 

Theorem 3.1. Thm.6.5,Thm.6.7] Consider PCR3BP with C = 3.03, fi = 
0.0009537. Then 

• there exist a homoclinic orbit to the L\ orbit (in Sun region). This orbit 
is close to the 3 : 2 resonance. 

• there exist a homoclinic orbit to the LZ, orbit (in exterior region). This 
orbit is close to the 1 : 2 resonance. 

These orbits are presented in Fig. ^ 

In this section we establish the existence of new homoclinic connections both 
in exterior and interior regions. The new homoclinic orbit in exterior region is 
close to the 2:3 resonance. As was mentioned in the Introduction this orbit has 
been found numerically in 0], see Fig. 5.4 the and the intersection stable and 
unstable manifolds of L* 2 at L — ^Ja ~ 1.26. The other new homoclinic orbit in 
interior region is close to the 5:3 resonance appears to be a new one. 



if P±(x )=xi 



then P±(R(xi)) = R(x ), 
then Pi i=F (J?(xi)) = R(x ). 
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Figure 1: 3:2 homoclinic orbit to L\ Lyapunov orbit (interior region) and 1:2 
homoclinic orbit to L\ Lyapunov orbit (exterior region). 
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Figure 2: 5:3 homoclinic orbit to L\ Lyapunov orbit (interior region) and 2:3 
homoclinic orbit to L\ Lyapunov orbit (exterior region). 
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3.1 The existence of the 2 : 3 homoclinic orbit in the 
exterior region. 

We define the following h-sets Gi = t(ci, Ui, Sj), for i = Q, . . . , 4, where 
c = (-1.12327231155833984,0), 

d = (1.093337837571255552,-0.02510094170679043584), 
c 2 = (1.047131544421841024,-0.001056187943513949696), 
c 3 = (1.08194053721089792,-2.521361165903333888- 10~ 5 ), 
c 4 = (1.04682616720451456,-9.169345277545603072- 10~ 7 ) 



and 



•so 


= (- 


1 • 10~ 8 ,4- 10~ 7 ), 


Uq = 


-R(s ), 


Si 


= (1 


• 10~ 7 ,21 • 10~ 8 ), 


Ul = 


-R( S1 )/10, 


S2 


= (- 


1 ■ 10- 7 ,35- 10" 8 ), 


U2 = 


-R{s 2 )/10, 


«3 


= (- 


1 • 10- 7 ,23- 10" 8 ), 


U 3 = 


-R(s 3 )/10, 


S4 


= (- 


1 • 10~ 7 ,35- 10~ 8 ), 


U4 = 


-i?(s 4 )/4. 



We assume, that Go,G2,G4 c 6+ and Gi,G 3 c 6-. With a computer assis- 
tance we proved the following 

Lemma 3.2. The maps 

Pi ;+ : GoUG 2 UG 4 ->e_, 
Pi,_ : Gi U G 3 6+ 

are u>eZ/ defined and continuous. Moreover, the following covering relations hold 

„ P l/2.+ _ -Pl/2 ; - _ -Pi/ 2,+ _ Pl/2,- _ Pl/2 ; + 2 

bO =?■ Ui =?■ 02 ==>■ 0-3 === ^ > ^4 =^ "2- 

Theorem 3.3. For PCR3BP with C = 3.03 and /z = 0.0009537 i/iere existe an 
or&ii homoclinic to L 2 close to the 2 : 3 resonance. 

Proof. From Lemma T3. 21 and |12l Lemma 5.6] it follows that 

„ Pi/ 2.+ „ Pl/2 - „ Pi/ 2 + n P i/2 - „ Pi/ 2 + jj-2 jP^ rr rr 
LrO =?* bl ( - r 2 =v- Lt3 => (_r4 => -f3 2 =?* "2 == ^ > -"2 

Note that the h-set Go is i?-symmetric by its definition. Therefore 

Hi = R(H 2 ) «=== R(H 2 ) «=== R{H%) P <^= i?(G 4 ) ^ + i?(G 3 ) 

i?(G 3 ) P <^= R(G 2 ) Pl ^= + Rid) Pl 4^ R(G ) = Go 

Since P_ is hyperbolic on \H 2 \ f |12l Lemma 5.5]) the assertion is a consequence 
of H Theorem 4]. □ 
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3.2 The existence of the 5 : 3 homoclinic orbit in the inte- 
rior region. 

As in the previous section we construct a chain of covering relations in order 
to prove the existence of homoclinic orbit to L\ orbit. We define h-sets Vt — 
t(ci,Ui, Si), for £ = 0, ... , 4, where 

c = (0.5217056203008400006,0), 

ci = (-0.5822638014577352639,-0.2793408708392046136), 

c 2 = (0.919204446847046941,0.004093829363524479834), 

c 3 = (0.9522506335647477061,0.0001333182992547130779), 

c 4 = (0.9208022956271231241,2.918364277340028028- 10~ 6 ) 

and 



so = 


(- 


1 ■ 10- 7 ,2 • 10~ 7 ), 


u = 


-R(s ), 


Si = 


(2 


• 10~ 8 ,4- 10~ 7 ), 


Ul = 


(3-10- 8 ,0), 


S2 = 


(- 


4- 10~ 7 ,102- 10~ 8 ), 


U2 = 


-R(s 2 )/5, 


S3 = 


(- 


1 • 10- 7 ,365 • 10~ 9 ), 


"3 = 


-R(s 3 )/io, 


S 4 = 


(" 


1 • 10- 7 , 25733011 • 10~ 14 ), 


Ui = 


-R(s 4 )/2. 



We assume, that Vo, V2, V4 c 6+ and V\, V3 C 0-. With a computer assistance 
we proved the following 

Lemma 3.4. The maps 

Pi t+ : V U V 2 UVi -► 6_, 

Pi_ : v 1 uv 3 ^e+ 

are well defined and continuous. Moreover, we have the following chain of cov- 
ering relations 

V V 1 Pl J^ V 2 V 3 P ^f- V, ^ H\. 

Theorem 3.5. For PCR3BP with C = 3.03 and /x = 0.0009537 there exists an 
orbit homoclinic to L\ close to the 5 : 3 resonance. 

Proof. Since the sets Hi and H\ are iJ-symmetric ^1 Lemma 5.6] and ^3 
Corollary 3.14] imply that 

H\ = R{Hl) R(H-l) =H 1 ^H 1 . 
After combining the above with Lemma 13.41 we obtain 

V P ^ + Vl P ^ V 2 ^ V 3 P ^ V, 2*> Hi ^ Hl ^ Hl . 
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Note that the h-set Vq is i?-symmetric by its definition. Therefore 



H X Hi R(V 4 ) P m R(V 3 ) R(V 2 ) P m W) ^ m) = V 

Since P+ is hyperbolic on \Hi\ f |121 Lemma 5.5]) the assertion is a consequence 
of H Theorem 4]. □ 



4 Symbolic dynamics on six symbols and reso- 
nance transitions. 

As a consequence of theorems proved in |12| and in the previous section we 
obtain the existence of symbolic dynamics on six symbols. Let L%, L 2 denote 
the Lyapunov orbits regions (see 0j) , S and / denote two parts of the Sun region 
corresponding to suitable vicinities of two homoclinic orbits to L* . Let X and E 
denote two parts the exterior region corresponding to suitable vicinities of two 
homoclinic orbits to L\ orbit. Schematically this situation is shown in Fig. [21 

In the symbolic dynamics on four symbols, i.e. {L 1 ,L 2 ,X, S} was es- 
tablished. The new homoclinic orbits allow us to include more symbols in it. 
We state this result more precisely. Let a, (3 £ {Li,L 2 , X, E, I, S} be such that 




Figure 3: The graph of symbolic dynamics on six symbols, 
there is an arrow from a to (3 on the graph presented in Fig. |31 We define the 
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function 



f/3,a = 



p+, 


if 


(a,0) 


= {L U L 1 ) 


P-, 


if 


[a, (3) 


= (L 2 ,L 2 ) 


P- ° (A/2.+ Pl/2,-) 4 P 1/2> + O P +! 


if 


(a,f3) 


= (Li,L 2 ) 


P+ o P 1/2) _ o (P 1/2 , + o P 1/2 ,_) 4 o P_, 


if 


[a, (3) 


= (i 2 ,ii) 


P+ ° (Pl/2,- ° Pl/2, + ) 2 ° Pl/2,- 


if 


[a, (3) 


= {S,L!) 


P+oP 1A _o(P lA+ oP lA „) 2 


if 


[a, (3) 


= (Li,S) 


P- o Pi/ 2 ,+ o (P 1/2 ,_ o P 1/2 , + ) 2 


if 


[a, (3) 


= (X,L 2 ) 


(Pi,_ o P 1/2 ,+) 2 o P 1/2j _ o P 2 


if 


[a, (3) 


= (L 2 ,X) 


P-°(Pl/2, + °Pl/2,-) 2 °Pl/2,+ , 


if 


[a, (3) 


= {E,L 2 ) 


Pl/2,- °(Pl/2,+ ° P/2,-) 2 Oil, 


if 


[a, (3) 


= (L 2 ,E) 


P-?°(Pl/2,-°Pl/2, + ) 2 , 


if 


[a, (3) 


= (I,Li) 


JPl/2,-°Pl/2,+ ) 2 oP-?, 


if 


[a, (3) 


= {LiJ) 



(4) 



For each symbol a 6 {Li, L 2 , X, E, I, S} we define the h-set Q a , where Qij = 
Hi, Ql 2 = H 2 , Qs = Eq, Qx — F , Qi = Vq, Qe = Go- 

Definition 4.1. The bi-infinite sequence (oij)jez is called admissible if for every 
i e Z there is an arrow from on to ctj+i on the graph presented in Fig. 

The finite sequence (ceo, cti, . . . , ct n ) is called admissible if for every i = 
0, 1, ...,n— 1 there is an arrow from on to cti+i on the graph presented in 
Fig.® 

Let r be the set of all admissible sequences (aj),gz € {Li, L 2 , X, E, I, S} z . 
Theorem 4.2. For every (cKj)j e z € T there exists a sequence (xt)i^z satisfying 
1. Xi E \Q ai \ for i 6 1, 
2- fai+i,an(Zi) = f or i G Z - 

Moreover, we have 

periodic orbits: if the sequence (o^igz is periodic with the principal period k 
then the trajectory may be chosen so that Xk — Xo, hence its trajectory 

is periodic 

homo- and heteroclinic orbits: if the sequence (c»!j)j£z is such that = Lj_ 
for k < fc_ and = L i+ for k > fc +; where € {1, 2} i/ien 



lim a^fc 

& — >oo 



»+' 



lim x/c = L*_ , 

/c — > — oo 

Proof. The same as ^1 Theorem 7.1]. □ 
4.1 Resonance transitions. 

Theorem 14.21 implies the possibility for a comet to move between various res- 
onances. If we interpret staying close to L* or L 2 periodic orbits as the 1 : 1 
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7 (a) e N le 



Figure 4: An h-set N and a horizontal curve 7 in TV. 



resonance, then Theorem 14.21 says that the comet can travel between exterior 
and Sun regions in both directions and can move between 5:3,3:2,1:2,2:3 
and 1 : 1 resonances in an arbitrary order. 

5 Symmetric periodic and homoclinic orbits. 

In Section |2~2*1 the symmetry property of PCR3BP and the associated Poincare 
maps are described. In this section we give the proof of the existence of an 
infinite number of symmetric periodic and homoclinic orbits. 

Definition 5.1. Let I 3 t — > u{t) G M 4 be a solution of PCR3BP, where I is 
the maximal interval of the existence of the solution. An orbit t — > u(t) is called 
i?-symmetric iff 

Image(u) = {u(t) | t G /} = {R(u(t)) | t G 1} = iZ(Image(«)). 

In this section we apply the method for finding symmetric periodic, homo 
and heteroclinic orbits first introduced in |1(J1 ITT] for the planar case and later 
developed in |13j in multidimensional situation. We recall here the basic defini- 
tions. 

Definition 5.2. Let N be a h-set with one unstable and one stable direction 
and let 7 : [a, b] —> R 2 be a continuous curve. We say that 7 is a horizontal 
curve in N if the following conditions hold: 

1. 7 ((o,6))cint(|JV|) 

2. either 7(a) € N le and j(b)eN re , 
or 7(a) G N re and 7(6) G N le . 

The geometry of this concept is shown in Fig. 

Definition 5.3. Let N be a h-set with one unstable and one stable direction 
and let 7 : [a, b] — > M 2 be a continuous curve. We say that 7 is a vertical curve 
in N if the following conditions hold: 
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1. 7 ((o,6))cint(|JV|) 

2. either 7(a) G N te and 7(6) G N be , 
or 7(a) G N te and 7(6) G N be . 

The following theorem is a special case of ^| Thm.3]. 

Theorem 5.4. Assume Nq , N\ , . . . , Nk are h-sets with one unstable and one 
stable direction and 

N a 4^N x 4^---^N k . 

If 7 : [a, 6] — > R 2 is a horizontal curve in Nq and 7 : [a, b] — > R 2 is a vertical 
curve in Nk then there exists to G (a, b) such that 

(fm o ■ ■ ■ o f o 7 )(t ) G int(|iV m+1 |), (5) 

/or m = 0, . . . , fc — 1 and 

(/fc-i°---°/o°7)(*o)G7((o,5)). (6) 

Theorem 15.41 was first proven in JU] for a planar case and direct covering 
relations. The generalization to a higher dimension with one unstable direction 
and the direct (forward) covering is presented in The proof of a general 

situation (i.e. direct and backward covering in multidimensional case) requires 
more sophisticated techniques and is presented in |13j. 

5.1 Symmetric periodic orbits. 

In this section we will use Theorem 15.41 in order to prove the existence of an 
infinite number of geometrically different symmetric periodic orbits. 

Before we state the main result in this section we introduce some notation. 
Let (a, (3) £ {Li, L2, X, E, I, S} 2 be an admissible sequence of symbols. Let the 
maps fp !a be defined as in Q. 

Notation: By Q a <=^t> Qp we will denote the chain of covering relations asso- 
ciated with the sequence (a,/?), i.e. 

where fp. a = P/- ° . . . ° Pi and Vi, i = 1, . . . ,k — 1 are suitable h-sets. 

Definition 5.5. Let f : X — > X. By Fix(/) we will denote the set of fixed 
points of f , i.e. 

Fix(/) = {y G X I f(y) = y}. 

Theorem 5.6. Let <\> : R x R 4 -o^R 4 denotes the local flow induced by the 
PCR3BP with C = 3.03 and n = 0.0009537. Assume (a , «n) S 
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left edge N 



Figure 5: A symmetric h-set. The 7 curve is both horizontal and vertical in N. 

{S, I, X, E, Li, L2}™ 7 n > is admissible sequence of symbols. Then there exists 
a point xo G \Qa \ such that 

(/a m ,a m _i ••• f ai ,a )(xo) G |Qa m |i 

(Ca^-^OWelQaJ, (7) 

for m = 1, ... ,11, i.e., the trajectory of Xq is coded by the periodic sequence of 
symbols 

(a n ,...,ai,ao,ai,...,a n ). (8) 
Moreover, xq is periodic and its orbit is R-symmetric. 

Proof. From the definitions of the h-sets used in the proof of homo- and hete- 
roclinic chains it follows that the sets H 1 , , Vq , Gq , Eq , Fq are i?-symmetric. 
Therefore Fix(i?) may be parameterized both as a horizontal and as a vertical 
curve in each of these sets (see Fig - Let 7 : [a,b] — > |Q ao | be the horizontal 
curve in Q ao and 7 : [a, b] — > \Q a „\ be the vertical curve in Q an , such that 
7([a, &]) U 7([a,6]) C Fix(-R). Now, Theorem 15 . 41 applied to the sequence 

Q f <^° O f ^> O f <^ 2 ■ ■ ■ fa !X^X- 1 n 
implies that there exists a point xq — 7(^0) G \Qa \ H Fix(i?) such that 

(/a m ,a m _i ' ' ' fa u a )( x o) G \Qa m \, for m = l,...,n, 
(/a„,a„_i o ■ • ■ o f auao )(x ) G j{{a,b)) C Fix(i?). 

From the definition of fp^ (see Eq.Q) as a composition of suitable Poincare 
maps it follows that there exists T > such that 

<p(T, Xq) = (/„„,«„_! ' ' ' fa u a„)(x ) G Fix(i?). 

Since i? is the reversing symmetry of 4> we obtain 

cj>(T,x ) = R(<f>{T,x )) = cj>(-T,R{x )) = <t>(-T,x ) 
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which proves xq is periodic and its orbit is i?-symmetric. 

There remains to prove that the trajectory of xq is coded by the sequence 
||SJ|, i.e. J7J) is satisfied. We formulate this as a separate lemma. □ 

Lemma 5.7. Assume (ao,...,a n ) is an admissible sequence of symbols. If 
x E dom(f olntan _ 1 o • • • o f ai . ao ) then R(x) G dom(/ Q T i 1 _ l:Ctn o • • • o f~^ ai ) and 

for m — 1, . . . , n. Moreover, if x — R(x) then 
for m = 1, . . . , n. 

Proof. One observes that if (a, (3) is admissible, then (/?, a) is admissible, too. 
Moreover, from the definition of /jg iQ (Eq. Q) it follows that Rofp,a — fa}j°R- 
Let x £ dom(f akt0lh _ 1 o • • • o f ai . ao ). Then 

(R° fa m ,a m -i ° ' ' ' ° fa u a )(x) = 
(/a„!-!,a m ° -R% m -i,a m - 2 ° ' ' ' ° /ai.ao)^) = 

for m = 1, . . . , n. If in addition x = R{x) then x G dom(/~^_ i a o • • • o f~^ ai ) 
and 

□ 

Remark 5.8. Theorem \5. 61 imvlies that there exist infinitely many geometrically 
different symmetric periodic orbits. This follows immediately from the fact that 
there exists an infinite number of admissible chains satisfying the assumptions 
of Theorem \5.b\ 

5.2 Symmetric homoclinic orbits. 

In this section we apply Theorem 15.41 in order to prove the existence of in- 
finitely many geometrically different symmetric homoclinic orbits to L\ and 
Lyapunov orbits. 

The following theorem shows how to use the method of covering relations 
in order to prove the existence of symmetric homoclinic or heteroclinic orbits. 
Later we will apply it to Poincare maps for PCR3BP. 

Theorem 5.9. Let No, Ni, . . . , N^ be h-sets, such that 
Nq 4=^ Ni 4=^ ■ • • 4=^ Nk Nk 
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and let 7 : [a, 6] — * M 2 in N be a horizontal curve in N . If fk is hyperbolic 
(see Def. 1]) on Nk, then there exists a point xq G j((a,b)) such that 

(fm o ■ • • o / )(x ) S int(|iV m+ i|), for m = 0, 1, . . .,fc - 1, 
{fk ° fk-l o ■ ■ ■ o / )(zo) G int(|JV fc |), for n > 0. 

Moreover, 

lim (/fc /fc-i ' ' ' fo)(xo) = 

n — >oo 

where x* is a unique fixed point of fk in \Nk\- 

Proof. From Theorem 15.41 it follows that for every n > there exists t n £ [a, 6] 
such that 

(fm ' ' ' fo)(l(t n )) G int(|iV m+ i|), for m = 0, 1, . . . , k - 1, 
(/?°/fc-i°---°/o)(7(U)Gint(|iV fe |). 

Since 7([a, 6]) is compact we can find G [a, 6] such that 

(fm° ■ ■■ °fo)(l(t*)) G int(|iV m+ i|), for m = 0,1,..., fc- 1, 
(fk ° /fc-i ° • • • ° /o)(7(*.)) G int(|AT fe j), for n > 0. 

Since neither /(7(a)) ^ iVi nor f(j(b)) £ N\ wc get G (a,o). Now, fk is 
hyperbolic on iVfc. Therefore by Theorem 3 in [2], 

lim (fk fk-i ° ■ ■ ■ ° fo)(xo) = £*, 

n — >oo 

where x := j(t*). □ 

Now we can state the basic result in this section. 

Theorem 5.10. Assume (ctQ, Ofi, . . . , ct n ) is an admissible nonconstant chain 
of symbols {S*, /, X, E, L\, L 2 }, such that a n € {Li,L 2 }. Then there exists a 
symmetric homoclinic orbit associated with the sequence of symbols 

(. . . , a n , a n , ctn-i, . . . , cti, ao, ai, . . . , a n -i,a n , a n , ■ ■ ■)■ (9) 

Proof. Let 7 : [a,b] — > |<5 Qo | be a horizontal curve in Q ao such that 7([a, b]) C 
Fix(i?). From Lemma 5.5 in ^2] it follows that P + is hyperbolic on |iJi| = Ql ± 
and P_ is hyperbolic on \H 2 \ = Ql 2 - Since a„ G {Li,i 2 } Theorem 15.91 there 
exists xq G 7 ((a, b)) such that 

(fa m ,a m -i o-'-o f auao )(xo) G int(|Q Q J), for m = 1, . . . ,ra, 
(/«„,«„ o ••• o / Ql , Qo )(x ) G int(|<3 a J), for fc > 0, 

, lim (/«„,«„ /«„,«„-! ' ' ' /ai.aoX^o) = 
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where L = L\ or L = L\ is a unique fixed point in |Q a „|. Since xq — R(xq) 
Lemma 15 . 71 implies that 

(/ai_!,a m ° ' ' ' ° fa \ ai )( x o) G int(|Q„J), for m = l,...,n, 

(/a*a B ° /«,._„«„ ° " " " ° /ao'aj^o) € mt(|Q a J), for fc > 0, 

lim (/~„ fc Q „ ° /£_ 1)OB o ■ • ■ o f- ) ai )(x ) = R(L) = L 

which proves that the trajectory of xq is a symmetric homoclinic orbit coded 
by the sequence of symbols @- □ 

Remark 5.11. Theorem \5. 1 ill imvlies that there exist infinitely many symmetric 
homoclinic orbits which are geometrically different. This follows immediately 
from the fact that there exists an infinite number of admissible chains satisfying 
the assumptions of Theorem \5.1(A 

6 Technical data. 

The computer assisted proofs of Lemma EOl and Lemma l!01 will be not discussed 
here. All ideas involved in such proof were presented in |12|. The C++ sources 
containing the rigorous numerical proof of Lemmas from |12| . Lemma \\\.2\ and 
Lemma 13.41 is available at [H]. 

The program uses the interval arithmetic and set algebra package developed 
at Jagiellonian University by CAPD group pQ. 

7 Appendix. Reading resonances from the tra- 
jectory in rotating frame. 

We describe the heuristic approach, which allows to read the resonance from 
the inspection of the trajectory in the rotating coordinate frame. 

We assume that Jupiter and comet move in the heliocentric inertial frame 
in the counterclockwise direction and the distance comet-Sun has well visible 
maxima or minima along the trajectory. This means that an approximate ellipse 
on which the comet is moving has nonzero eccentricity. 

Let R denote the resonance. Let T be an approximate period of the comet 
in the heliocentric frame. Let us recall that the period of the Jupiter is equal 
to 1. Hence 

R=^ (io) 

Then in the heliocentric inertial frame the average angular velocity of the comet 
is and that of the Jupiter is equal to 2ir. 

For an approximate periodic trajectory of a comet in the rotating frame we 
introduce the following notation 
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9 is the number of full turns around the Sun during the whole period. 
This number is positive for trajectories in the interior region and negative 
in the exterior region. 

M - the number of maxima (or minima) of the distance between the Sun 
and the comet. 



Since the distance Sun-comet reaches the maximum (or minimum) only when the 
comet is at the aphelion (or perihelion), hence consecutive maxima (minima) 
occur with the period T. In the rotating frame the difference between the 
angular variables of the comet and Jupiter is equal to 2%6/M. Observe that 
this difference is the same in both reference frames, the inertial one and the 
rotating one. Hence 



2ir9 


/2tt 


~M~ 


\T~ 


6 


l-T 


M 






M — 9 


T = 






M 



Hence finally 

R = -^- a . (11) 

M -9 v ' 

Let us apply i|ll|) to Figures ^ and |2 For interior homoclinics we count the 
maxima and for exterior homoclinic we count the minima. We have 



the interior homoclinic orbit in Figure ^ 9 = 1, M = 3. Hence R = ~ = 
3:2. 

the interior homoclinic orbit in Figure 9 = 2, M = 5. Hence R = | = 
5:3. 

the exterior homoclinic orbit in Figure ^ @ = — 1 ; M = 1. Hence R = 
1 = 1:2. 

the exterior homoclinic orbit in Figure |2 = — 1, M = 2. Hence R = 
| = 2:3. 
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